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We propose a fast and non-destructive spectroscopic method for single molecular ions that im- 
plements quantum logic schemes between an atomic ion and the molecular ion of interest. Our 
proposal relies on a hybrid coherent manipulation of the two-ion system, using optical or magnetic 
forces depending on the types of molecular levels to be addressed (Zeeman, rotational, vibrational 
or electronic degrees of freedom). The method is especially suited for the non-destructive precision 
spectroscopy of single molecular ions, and sets a starting point for new hybrid quantum computation 
schemes that combine molecular and atomic ions, covering the measurement and entangling steps. 
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I. INTRODUCTION 

Recent advances in trapping and cooling of molecules 
and molecular ions are opening new perspectives in fields 
as varied as cold chemistry, quantum dynamics of com- 
plex systems, metrology and quantum computation [Ti- 
ll]. In this context, the sympathetic cooling of molecular 
ions represents a particularly attractive approach which 
enables to perform chemical and spectroscopic experi- 
ments on the single-particle level [5]. With the recent 
breakthrough in the production of cold molecular ions 
in well-defined internal states [6H2], fully coherent ex- 
periments with single molecular ions now come within 
reach. Precise knowledge and control of the internal 
quantum state of molecular ions is relevant for quantum 
controlled chemistry [3 [TO, frequency metrology of single 
ions [TTJ [12] and prospective applications in quantum- 
information processing, where a lot of progress has al- 
ready been made with trapped atomic ions [T3j [14] . How- 
ever, readout procedures of the internal quantum state 
have thus far relied on methods such as laser-induced 
charge transfer [8j [9] or photodissociation [6] [7] because 
of the lack of closed cycling transitions in molecular sys- 
tems. These techniques are inherently destructive and 
not suitable if repeated measurements on a single molec- 
ular ion are required, for instance in quantum computa- 
tion [14] or frequency metrology [TTj . 

In this work we propose a fast, efficient and accurate 
molecular spectroscopy scheme based on quantum logic 
and coherent control theory, able to address a variety of 
internal degrees of freedom of the molecules. Our pro- 
posal relies on a hybrid manipulation of the ion of in- 
terest and an atomic ion, such as 40 Ca + or 9 Be + , that 
acts as a probe. Optical or magnetic forces are used on 
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the ions at will, depending on the precise internal states 
which are involved in the quantum operations. They are 
arranged in optimized gates that take a time from 50 /is 
up to 1 ms, are insensitive to the temperature of the 
ion crystal and have a number of relevant applications. 
The most immediate one is the determination of mag- 
netic moments or Zeeman shifts in atoms and molecules, 
such as 14 N2~ or 16 Oj, for precision-spectroscopic pur- 
poses. This pushes the field along the line of previous 
experiments with atomic ions [11] , incorporating two im- 
portant features: a greater versatility (that allows ad- 
dressing a broader variety of ions with a complex inter- 
nal structure), and simpler experimental requirements, 
as cooling the ions to the motional ground state is not 
necessary. Another application of this work is to replace 
the destructive measurement techniques in current exper- 
iments with hybrid atom- molecule systems [HJ [9] . More- 
over, the proposed framework represents the corner stone 
of a hybrid quantum computation scheme which for the 
first time combines molecular and atomic ions, covering 
the measurement and entangling steps. Finally, in the 
broader context of molecular physics, the most impor- 
tant aspect of the scheme is that it represents an entirely 
new approach for molecular spectroscopy which relies on 
the manipulation of quantum phases and is tailored to- 
wards the interrogation of single trapped particles. 

Our manuscript is structured as follows. We start in 
Sect, [n] by reformulating the concept of quantum logic 
spectroscopy [11 , showing how by means of control-phase 
gates and Ramsey interferometry, a single ion may probe 
the state population or other observables of any other 
ion in a Coulomb crystal [5]. We then show how to im- 
plement those control-phase gates using state-dependent 
forces [T5HT7] in a crystal with different types of atomic 
and molecular ions, thus introducing a quantum proto- 
col for molecular ions. With these tools, in Sect. |III| we 
derive three slightly different variants of this protocol, 
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and test them with accurate calculations for the case of 
14 N2~ ions. In the first version, the probe and the target 
ions are subject to pulsed optical forces generated by AC 
Stark shifts. This fast protocol is targeted at the deter- 
mination of electronic, rotational and vibrational states 
of the target molecular ion. In the second version, the 
ions are subject to oscillating magnetic fields, allowing 
the determination of Zeeman states. The final version 
relies on femtosecond laser pulses acting on the control 
ion, leading to an order of magnitude greater sensitiv- 
ity of the detection scheme. We draw our conclusions in 
Sect. El 



II. QUANTUM LOGIC SPECTROSCOPY 
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Roughly speaking, spectroscopy studies the reaction 
of a physical system under periodic drivings, forces, or 
radiation of various frequencies. The most common ap- 
proach is to watch the back-action of the system on the 
radiation that drives it, measuring the absorption, the 
emission or the phase in those fields. The field of quan- 
tum logic spectroscopy (QLS), introduced in Ref. [TT] . 
advocates a very accurate study of the changes suffered 
by the driven system, using quantum gates to enhance 
the precision of those studies. In this section we extend 
these ideas to the field of molecular spectroscopy. To 
this end, we have to overcome the two main limitations 
of the original QLS protocol, which are the temperature 
requirements and its speed. 



A. Phase sensitive QLS protocol 

The original experiment for QLS [11 employed two 
atomic ions that allow optical manipulation, 9 Be + and 
27 A1 + . By means of sideband transitions, the state of 
the target ion, 27 A1 + , which cannot be directly mea- 
sured with spectroscopic accuracy because it lacks cy- 
cling transitions, was mapped onto the control ion, 9 Be + , 
which was later on measured using the accurate tech- 
nique of electron shelving. Hence, in this pioneering 
work, quantum logic addressed the lack of cycling tran- 
sitions in 27 A1 + , at the price of imposing accurate co- 
herent operations (sidebands) and ground state cooling 
of the motional state of the ions in the trap. Both re- 
quirements are experimentally demanding and have been 
addressed in later experiments with atomic ions (see 
e.g. [HJ [19]). However, they have not yet been achieved 
for sympathetically-cooled molecular species. 

One can dramatically enhance the versatility of the 
QLS protocol by using geometric gates based on state- 
dependent forces [T5HT7] . These are quantum gates that 
are insensitive to the motional state of the ions and 
may thus be applied to sympathetically cooled molecular 
ions. The underlying physical principle of these gates, ex- 
plained in Sect. |IIBl is that when a chain of ions is shaken 
by some external forces, the quantum state acquires a 



Figure 1. State detection protocol, (a) The probe ion experi- 
ences a series of local gates combined with a two-qubit control- 
phase gate in which the control is a property, O, of the target, 
(b) Oscillations in the state of the probe ion as a function of 
the acquired phase, <E>o oc tO. We plot the mean excitation 
probability with the error after 10 repetitions (£ = 0). 



phase that depends on those forces. Moreover, when 
the forces experienced by the ions are "state-dependent" , 
that is, when they have different magnitudes or signs 
depending on the internal state of the ions, then the 
quantum phase carries precise information on the state 
of those atomic or molecular ions. This information is 
the one we use for the QLS algorithm. 

More precisely, assume that we have two ions, a control 
or logic ion, labeled C, and a target or spectroscopy ion, 
labeled T, confined in the same ion trap. The C ion is 
Doppler laser cooled and sympathetically cools the tar- 
get ion to form a two-ion Coulomb crystal. Both ions will 
be subject to independent external forces, different from 
the laser-cooling process. The ion C will be a qubit, and 
the force acting on it will depend on its internal state, 
fc(t)o-Q, while the force on T will depend on some prop- 
erty of this ion, /t(^)0, such as a magnetic moment or a 
quantum number. Here, cfq is the Pauli ^-matrix for the 
internal state of the C ion. As explained in Sect. II B[ 
the quantum state of the system will acquire a total non- 
trivial phase of the form 



$o = crgO^cT, 



(1) 



where (/>ct is a function of the normal modes of the two 
ions, u com and u; str [20] , and their combined drivings. 
This phase is robust, it does not depend on the motional 
state of the ions at the beginning of the gate operation 
and thus is independent of temperature. Based on this 
we design a quantum protocol which interrogates very 
accurately the state- dependent phase <&o-> obtaining the 
value of the observable O (cf. Fig. [l^i) 

1. Prepare ion C in state |0); 

2. apply a Hadamard gate on C, H = exp(— icr^7r/2), 
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3. optionally, apply a reference phase on the control 
ion, exp^cr^); 

4. apply the state dependent forces /c,t(£) on the C 
and T ions, ensuring that the initial motional state 
of the ions is restored (see Sect. IIIB); 



5. undo the Hadamard gate, iff; and 

6. measure the state of the control ion. 

At the end of this process the probe and the target ions 
will get entangled and the excited state population of the 
control ion will oscillate as (cf. Fig. [TJd) 



P t = sin 2 (£O0 C T+O- 



(2) 



This allows, by repeated application of the protocol, to 
determine the value of O with very high precision [2T] , 

We have so far presented the protocol in a very general 
form, without clarifying the origin and the time depen- 
dence of the forces that act on the control and target ion. 
In the following section we will elaborate on the physics 
behind the phase gate (TJ, developing the formalism to 
compute and optimize the control and target forces. Us- 
ing this we will then work out the actual implementation 
of the protocol to detect molecular electronic, rovibra- 
tional or Zeeman states, using optical or magnetic fields. 



B. Geometric phases 

As pointed out in Ref. [17] , quantum harmonic oscilla- 
tors possess a very simple geometric structure: a continu- 
ous displacement of the oscillator through phase space in- 
duces a global phase on the wavefunction, and this phase 
only depends on the initial and final points of the path, 
i.e., it is geometric in nature. Following the framework in 
Ref. [15 , one can show that when a harmonic oscillator 
experiences a force hf(t) for a time T, its quantum state 
acquires a phase 

<j>[u, f]=al I [ sin[cj(ri - T 2 )}f(T 1 )f(T2)dT 2 dT 1 , 
Jo Jo 

(3) 

where uj and ao are the harmonic oscillator frequency and 
length, respectively. This phase is proportional to the 
area in phase space covered by a coherent state subject 
to this force, whose trajectory is 

z(t; u, f) = z(0) -i^f f(r)e^dT. (4) 

This geometric phase is resilient to temperature, for the 
area does not depend on the initial state of the oscillator, 
z(0). Moreover, it is also insensitive to random errors in 
/(£), which only appear as second-order corrections. 

In our protocol we will act with two independent forces, 
fc and /t, on two ions. In this case we have not one but 
two harmonic oscillators, associated to the center of mass 



(com) and stretch mode (str) of the ion crystal. Since the 
target and control ions are different, they will in general 
have different masses and also experience different con- 
finement potentials. This complicates the final expres- 
sions, but we can still write the effective Hamiltonian 



H 



sG {com, str} 



hw s a\a a 



■hF s (t) — (a s 



4)} 



(5) 



where F com and F str are linear combinations of cr^fcit) 
and 0/t(£) (cf. Appendix [A]). The total phase acquired 
by the combined system may be written 



$ = </>[aWi, ^com] + </>[^str, F 8tr \. 



(6) 



In this expression there will be many contributions, but 
the only one that will influence the signal of our protocol 
is the one proportional to cfq x O, for it is the only en- 
tangling phase in this process. This nonlocal phase may 
be rewritten as the product ([I]), with some integral (/>ct 
that depends on the actual driving, but which may be 
computed analytically or numerically. 



III. IMPLEMENTATION OF THE QLS 
PROTOCOL 

In the following sections we will present three different 
versions of the protocol that can use different types of 
forces and drivings, and how they can be used to probe 
different properties of the target ion. In order to derive 
analytical estimates, we will work with Gaussian pulses 
having a duration T and driven at frequency v 



/c,TW = /8, T e- (2 * /T)2 cos(^). 



(7) 



This ansatz does not represent any loss of generality, as 
the quantum gates may be further optimized choosing 
more sophisticated dependencies. Nevertheless, experi- 
ence shows that the limits set by the previous pulses are 
very close to those achievable by more general forces [16] . 

For this time dependence of the driving, and a suffi- 
ciently long pulse, T > 57r/min{|^ — o; com |, \ v — u; str |} so 
that all motional degrees of freedom are restored to their 
initial state, the two-ion state accumulates a total phase 
as given by Eq. with (cf. Appendix |A|) 



1 /^ /g/°a 2 T 
^ CT = 4V2^^ (W) - 

Here we introduced the enhancement function 



E(i/,u;,/i) := 



UJ 



UJ 



(8) 



(9) 



where \i = mx/^c is the mass ratio of the two ions, 
a 2 = H/[(mc + %)4 and uj is the frequency of the 
control ion trap. Note that our formalism allows the 
study of arbitrary drivings, ^, which can be advantageous 
in some setups, as we discuss below. 
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Figure 2. (color online) Scheme for imparting an optical force 
on an ion. (a,c) A single off-resonant laser beam is focused on 
the particle with a radial displacement, so as to achieve the 
maximum intensity gradient. (b,d) Two counter-propagating 
laser beams create an optical lattice and the atom sits on the 
point of maximum slope. 



A. Optical kicking protocol to address electronic 
and rovibrational states 



The first version of the protocol aims at distinguishing 
electronic, vibrational and rotational states of the molec- 
ular ion. For these degrees of freedom it is both advan- 
tageous and very efficient to rely on optical "kicks" as 
the state-dependent driving mechanism behind the gate. 
However, due to the inherent complexity of the molecu- 
lar energy- level structure, we must seek gentle ways to 
drive the charged particles without accidentally exciting 
any undesired internal transitions. For this purpose, we 
will employ the dipole forces induced by two lasers that 
are detuned from an electronic transition of the atomic 
and the molecular ion, respectively. Each ion will sit on a 
spot of maximum slope of intensity, experiencing an AC 
Stark shift force 



/c,tW 



n(t) 2 



1 

T 



(10) 



where A is the laser detuning, Q is the Rabi frequency 
and Vt 2 is proportional to the light intensity [22]. The 
force depends on some characteristic length over which 
the light intensity varies, £, which ranges from half a 
wavelength for an optical lattice, up to £ ~ 5 — 100/im 
for a focused laser beam (see Fig. [2|. 

We will work in the pulsed regime, with v = and 
T <C 27t/cj, applying four kicks on the ions. The kicks 
will be spread in time according to 

{(Afc, -ti), (Afc, -t 2 ), (-Afc, * 2 ), (—Afc, *i)}, (11) 

with (ti,t 2 ) = (0.920, 0.080)2tt/cj, see Fig. ^ This se- 
quence, which takes two periods of a trap, is similar to 
the one developed in Ref. [16 for femtosecond resonant 
laser pulses, but the forcing mechanism in that reference 
can not be implemented with molecules because spec- 
troscopic addressing is very hard. We emphasize that 
we rely instead on ultra-fast off-resonant kicks. Interest- 
ingly, as we show below, working with off-resonant kicks 
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Figure 3. (color online) Effect of the four pulses (11) on a 
coherent state of a harmonic oscillator. The red dashed lines 
correspond to the unforced motion while the solid blue lines 
are subject to the kicks, (top) Phase-space plot: the trajec- 
tory starts at the point indicated by the small black circle. 
The arrows indicate the effect of each successive pulse, (bot- 
tom) Position, Q, and momentum, P, as a function of time, 
t. Note how both variables return to their initial orbit after 
the fourth kick. 



does not result in particularly demanding field intensi- 
ties. 

The momentum kick can be approximately related to 
the desired phase ([I]) as follows 



Afc ~ TJf°a c fta T ~ CT- 



(12) 



Requiring that the photon scattering probability remains 
below the small value e — TQ?T / A 2 <C 1, we can extract 
the product of the average Rabi frequency, and the 
kick duration, T, as a function of the detuning, A, the 
spontaneous emission rate of the ions, Tc,t? and the error 
tolerance, e. From here we obtain the relation 



CT 



A a c ,T 

£X r x — 



(13) 



which states that we need low trap frequencies, steep 
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Figure 4. Optical kicking protocol acting on a 40 Ca + and a 
14 Nj ion. (a) Required detunings vs. trap frequency, for a 
scattering error e — 1% and beam waist £ = 5/xm (solid), 
and e = 0.1% in an optical lattice, £ — A/2 = 397 nm 
(dashed), (b) Gate pulse time (solid) and total gate du- 
ration including pauses (dash-dot), for e — = 5/xm 
and 1 mW laser power. For simplicity, we assumed the 
same detuning for the lasers acting on the 40 Ca + 
2 P 1/2 , 397 nm) and 14 N+ (|X 2 S+, v" = 0, J" = 1/2, Fi) 
|£ 2 E+, v = 0, J' = 3/2, Fi), 391 nm) transitions. 



light gradients and large detunings in order to acquire a 
large phase shift. 

We have studied numerically the conditions to per- 
form these quantum gates with a variety of atomic and 
molecular species. As an example, in Fig. we plot 
the required detunings for making an accurate quan- 
tum gate (e ~ 0.1%, dashed) and for running the spec- 
troscopy protocol (e = 1%, solid) using a setup with 
40 Ca + as the control ion and as the target ion. 

For the atomic ion we chose the cycling transition that 
is excited in electron-shelving measurements, 2 S1/2 ~^ 
2 P1/2, for which existing lasers may be conveniently 
reused. For the molecular ion we choose the strong 
electronic transition |X 2 £+, v" = 0, J" = 1/2, Fx) 
\B 2 Z+,v' = 0, J' = 3/2, Fi), which would allow us to 
detect states that are already produced in the labora- 
tory [8 . As expected, if we need greater accuracy in 
the two-qubit gate without exceedingly large laser power 
(e ~ 0.1%), it is best to adopt an optical lattice configu- 
ration with the ion sitting on the slope, cf. Fig. [2}},d. If 
we are only interested on spectroscopy and state deter- 
mination, we may increase the error tolerance to e ~ 1% 
and place the ions at half the waist from the center of 
a focused laser beam, cf. Fig. [2^i,c. In both cases the 
duration of the pulse is much shorter than the trap fre- 
quency, and gate times of microseconds are feasible using 
laser powers of milliwatts. 
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Figure 5. (color online) Zeeman energy shifts of the hyper- 
fine levels of the ground rovibronic state of 14 (left) and 
their derivatives, \i — (dE/dB)/K (right), which are needed 
to estimate the state dependent force (14). 



B. Magnetic forces for Zeeman states 

In order to discriminate Zeeman states within a molec- 
ular rovibrational level, we develop a second version of 
the protocol, based on the use of magnetic forces, which 
result from the coupling of a magnetic field gradient with 
the effective magnetic moment /i of the particle [23j [24] , 



hit) 



dB(t) n 
dx h 



(14) 



In this setup spontaneous emission is not a concern, but 
the forces are going to be typically much weaker. This 
may be compensated by driving the magnetic field close 
to resonance with one vibrational mode of the ion crystal 
as realized in Ref. [25]: for a driven force [Eq. ([7])] with 
v = (1 + T])u 

com close to the center of mass mode, the 
phase is enhanced, $(77) ~ ^(0) / 77, at the expense of a 
longer gate, 1 / 77 times longer. 

As a relevant, realistic example, we study the Zeeman 
splitting of the hyperfine structure of the 14 Nj molec- 
ular ion in the rovibrational ground state of the elec- 
tronic ground state potential X 2 S+ [3|9]. The molecu- 
lar Hamiltonian in the presence of an external magnetic 
field reads 



(15) 



H — H rot + H sr + i^hfs + H eq Q + Hz , 



where H rot is the rotational Hamiltonian, H sr is the spin- 
rotation interaction, i^hfs is the magnetic hyperfine in- 
teraction, neglecting the coupling between the external 
magnetic field and the nuclear spin, H eq Q is the nuclear 
electric quadrupole interaction, and finally Hz is the Zee- 
man interaction. We have calculated the lowest Zeeman 
levels and their respective magnetic moments, shown in 
Fig. [5] (see Appendix [B] for details). It is remarkable that 
almost all states are distinguishable by their magnetic 
moments, not requiring more than 10% accuracy in its 
determination. 

Using the magnetic moments one can easily compute 
the field intensities required to perform interferometry 
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Figure 6. (top) Gradient and (bottom) time required to 
achieve maximum entanglement between 40 Ca + and 14 N J 
ions, using an oscillating magnetic field gradient dB/dx = 
exp[— (2t/T) 2 ]B / cos(z/£) in a harmonic trap, cj, for drivings 
v/oj = 0, 1.1 and 1.01 (dash-dotted, dashed, solid). 



using a 40 Ca + and a 14 N+ ion, as a function of the trap 
frequency. As shown in Fig. |6j a trap with uj = 2tt x 
574 kHz requires a magnetic field gradient of 10 T/m for 
a 250/is gate using a driving v = I.OIcj, values which do 
not seem disparate [25, 29 j. However, note that for the 
same duration there exists another configuration for v = 
which produces the same gate at a lower value of B' = 
dB/dx, using a smaller trap frequency. More precisely, if 
we can make a gate with field gradient B r using v = (1 + 
r/)w com in a time T v ~ bn/rjuu, we can perform the same 
gate using a smaller trap frequency uj u =o ex rf^uj v and 
a significantly shorter time T u=0 ex r] 1/3 T u . This setup 
may be advantageous for larger traps, where achieving 
a large B' may be difficult and would otherwise require 
long gate times; using in this case v = would reduce 
the gate time, thus avoiding heating problems [301 131]. 



C. Opto- magnetic protocol 

Both of the preceding methods may hit an important 
problem with certain molecules which due to their mass 
or their internal structure respond more weakly to the 
external optical and magnetic forces. For such cases, we 



propose a third protocol that combines the application 
of ultrashort laser pulses on the control ion with arbi- 
trary optical or magnetic forces on the target. Driv- 
ing the control ion with a sequence of N femtosecond 
7r-pulses [16] we will be able to inject an almost arbitrary 
amount of momentum Akc ~ N x 2tt/X to the system. 
The enhanced pushing capability of the pulse train com- 
pensates the smaller kicks, Afcx? suffered by the target 
particle, as the total phase is proportional to the prod- 
uct, (J)ct ~ A&cA&t. This simple idea may be used in an 
ultrafast combined opto-magnetic protocol to lower the 
magnetic field requirements by various orders of magni- 
tude 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we have introduced a spectroscopy pro- 
tocol based on quantum logic that is applicable to a wide 
range of ions without accessible cycling transitions, such 
as molecular ions. We have presented numerical results, 
based on original molecular structure calculations, show- 
ing that the requirements to implement this protocol with 
molecular ions of current interest are within experimental 
reach. This protocol, which combines magnetic and op- 
tical forces, enables the realization of measurements and 
entanglement on hybrid atom-molecule ion systems with- 
out the need for cooling to the trap vibrational ground 
state, a significant improvement over previous methods. 
This constitutes a new approach to high-resolution spec- 
troscopy of cold molecular ions, and opens the door to- 
wards novel hybrid quantum computation schemes. 



ACKNOWLEDGMENTS 

This work is supported by Spanish MICINN Projects 
FIS2009-10061, FIS2010-22064-C02-02 and CTQ2007- 
62898-BQU, CAM research consortium QUITEMAD 
S2009-ESP-1594, the Swiss National Science Foundation 
through the National Centre of Competence in Research 
"Quantum Science and Technology" , ESF COST Action 
IOTA (MP1001), a FP7 Marie Curie fellowship (IEF- 
2009-251913 MOLOPTLAT), and a JAE CSIC Fellow- 
ship. 



Appendix A: Calculation of the total phase acquired 

We sketch here the derivation of the phase acquired by 
the two- ion state due to the forces applied on the ions, 
Eq. We start by recalling the total phase acquired 
by a state of a harmonic oscillator of frequency uo driven 
by a force f(t) for a time T [T5] : 

<)) =l^l m j\tj\Te i ^ t -^f{t)f{T). 
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For a time dependence as in Eq. ([7]), and assuming a suf- 
ficiently long pulse so that the motional state is restored 
to its initial value, T > 57r/min|^ — cj|, this phase takes 
the form 



fa 

h 



ujT 



8(w 2 - I/ 2 ) ' 



(Al) 



with a — yjhj (muj) the harmonic oscillator length. 

For the case of a two-ion system, we have to take into 
account that a force applied on one of the ions will af- 
fect the two common modes, center-of-mass (com) and 
stretch (str). Applying forces F\ and F 2 on the two ions, 
the forces on the com and str modes are 



^com — F\ + F 2 , 



mi + rn 2 



(m 2 Fi - miF 2 ) . 



The p hase acquired is then given by Eq. (|6| . Using the 
result (Al), we obtain that the only contribution to 



that depends on F\F 2 is 
T 

$ F X F 2 



2 8(mi + m 2 )h W om ~ y2 ^str ~ 1 
which is readily rewritten as Eq. ([8| with T = F\^ 2 /h. 



Appendix B: Calculation of 14 Nj hyperfine levels 



The complete molecular Hamiltonian for a 2 £ molecule 
such as Nj, including the hyperfine structure and the 
corresponding splittings under a magnetic field £?, in- 
volves several angular momenta: the rotational angular 
momentum iV, the (total) electronic spin S, and the to- 
tal nuclear spin I. These momenta are coupled within a 
Hund's case (b) scheme, i.e., J = N + S and F = J + I. 
The terms in Eq. ( 15 ) are given by [26-28 



• The rotation of the molecule including the centrifu- 
gal distorsion, 



rot 



B P N Z - D P N 



(Bl) 



• The fine structure spin-rotation interaction, 

H sr = 7sr iV • S. (B2) 



• The magnetic hyperfine Hamiltonian that couples 
the nuclear spin with the electronic spin and molec- 
ular rotation, 



#hf s = H 



is 



H 



IN- 



(B3) 



The first contribution can be written as a sum of 
the Fermi contact interaction, = bpI-S, and the 
dipolar term, H dip = -VlOT^(I) -T^(S,C 2 ), 
in spherical tensor notation (see Ref. [28] for de- 
tails). The second part, ii/jjv, is similar to the 
spin-rotation term but considering the nuclear spin 
instead, Hin = cjl • N. 

The nuclear electric quadrupole interaction, which 
can be expressed using spherical tensor notation as 



(B4) 



with e the electron charge. 

• Finally, for the Zeeman interaction we include only 
the coupling of the electronic spin with the mag- 
netic field 



Hz g">BBS z 



(B5) 



where g is the electron gyromagnetic factor (taken 
g = 2), B is the magnetic field, hb is the Bohr mag- 
neton, and S z is the component of the electronic 
spin along the quantization axis that is taken in 
the direction of the external field. Any other mag- 
netic contributions will be neglected as they are of 
much less importance. 

We work within the Hund's case (b) and select the 
basis \N SJIFMp) to represent the Hamiltonian. The 
nuclear spin of 14 N equals I\ = 1, and the calculations 
have been carried out for the symmetry block which in- 
volves even values of both the rotational (N = 0,2,4,...) 
and the total nuclear spin (I = 0,2). For zero magnetic 
field, our results are in complete agreement with those 
appearing in [26 . We have also checked that the calcula- 
tions reported here are converged just including N = 0, 2 
in the Hamiltonian matrix. 

With the values for the parameters of the molecular 
Hamiltonian reported in Table [l| the matrix elements 
read explicitly as follows [26] : 



• The rotational term: 

{N' S r J' I' F' M' F \H vot \N S J I F Mp) = ^nn'^ss'^jj'^ii'^ff'^MfM' 

(B e N(N +1) - D e [N(N +1)} 2 ) . (B6) 

• The fine structure, spin-rotation term: 

(N'S'J'I'F'M' F \H sr \NSJIFM F ) = 5 NN ,6ss'Sjj,6 n ,5 FF ,6 MFM ^ [J(J + 1) - N(N + 1) - S(S + 1)} . (B7) 
Here, the coupling constant is defined by 7 sr = 7 + JnN(N + 1). 
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Table I. Values of the molecular Hamiltonian parameters used in the calculations. B e and D e are taken from Ref. [32] and are 
here given in cm -1 , while the other parameters are taken from [26] and given in MHz. 



Parameter 


B e D e 7 


7iv 


b F Cdip cj eqQ 


Value 


1.9223897 5.9758 x 10" 6 276.92253 


-3.9790 x 10" 4 


100.6040 28.1946 0.01132 0.7079 



• The magnetic hyperfine structure terms: 

{N'S'J'I'F'M' F \H ¥ \NSJIFM F ) = 8 NN ,5ss>hr5FF>5 MF M>M- l ) F+I+J+J ' +N+S+1 



x + 1)(2J + l)y/S(S + 1)(2S + lW(2J' + 1)(2J + 1) 

{N'S'J'I'F'M' F \H^\NSJIFM F )=5ss'hr5 FFI 5 MF M' F c A ^{-l) F+I+J+N ' +1 



{■ 


J' 






J' 






I 






s 





y/30 1(1 + 1)(2J + 1)S(S + 1)(25 + 1) V(2J' + 1)(2J + 1)(2N + 1)(2N' + 1) 





j' 






' N' 


N 


2 




F * 




S 










\\ 


S 


1 > 


I 












J' 


J 


1 



N f 2 N 




{N'S , J , I'F , M' F \H IN \NSJIFM F ) = 5 nn ,5ss>5ii>5ff>5m f m>M- 1 ) 



F+I+J' + J+N+S+l 



x + 1)(2/ + l)y/N(N + 1)(27V + l)y/(2J' + 1)(2 J + 1) 





1 






i 












s 





(B8) 



(B9) 



(B10) 



with 



>, 3-j, 6-j, and 9-j symbols, respectively. All matrix elements couple func- 



tions with AJ = 0, ±1. In addition, the dipolar interaction, i^dip? mixes states with A7V = 0, ±2. The values 
for ftp? c^p = t + ttfN(N + 1), and cj are given in Table |lj We have neglected the centrifugal distortion term 
tjsf because of its very small value. 

• The nuclear quadrupole term, that allows for a number of off-diagonal elements (A J = 0, ±1, ±2; AI = 0, ±2; 
AN = 0,±2): 



{N'S'J'I'F'M' F \H eqQ \NSJIFM F ) = 5 ss >Sff>Sm f 



eqQ + (-if 1)F+r+2J 

'2 2 1 ' 



(-1) 



2/i+S+2AT' 



x y/(2I + 1)(2J' + 1)(2J + 1)(2J' + 1)(2JV + 1)(2N' + 1) 



r 2 i 

J F J' 



h 2 h 
I h I 1 



N' 2 N 
J S J' 



N' 2 N 




h 2 h 
-h h 



(Bll) 



where /i is the nuclear spin of one of the nuclei composing the (homo nuclear) molecule, and the value of eqQ is 
also from Ref. [26] . 

• Finally the Zeeman term [28] : 

(N f , S', J', I', F', M' F \H Z \N, S, J, I, F, M F ) = 5 NN ,5 ssl 5 n ,g^ B B{-l) F ' ~ M '- ^J'+n+s+f+i 



y/(2F' + 1)(2F + 1)V(2J' + 1)(2J + l)y/S(S + 1)(25 + 1) 

( F' 1 F 
\ -M' F M F 



{' 


J' 






s 






F 


:) 




j 





(B12) 
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